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■ ABSTRACT 

Aims. Horizontal branch (HB) stars in globular clusters offer us a probe of the mass loss mechanisms taking place in red giants. For 
. M3 (NGC 5272), in particular, different shapes for the HB mass distribution have been suggested in the literature, including Gaussian 

' and sharply bimodal alternatives. In the present paper, we study the mass distribution of HB stars in M3 by comparing evolutionary 

tracks for a suitable chemical composition and photometric observations. 

Methods. Our approach is thus of a semi-empirical nature, describing as it does the mass distribution that is favored from the stand- 

■ point of canonical stellar evolutionary predictions for the distribution of stars across the color-magnitude diagram. More specifically, 
we locate, for each individual HB star in M3, the evolutionary track whose distance from the star's observed color and magnitude is 
a minimum. We carry out artificial tests that reveal that our method would be able to detect a bimodal mass distribution resembling 
the one previously suggested in the literature, if present. We also study the impact of different procedures for taking into account the 
evolutionary speed, and conclude that they have but a small effect upon the inferred mass distribution 
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1 . Introduction mal deviates were assumed [after Rood & Crocker, who advo- 

cated a normal distribution for the masses of horizontal branch 
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^ ' Results. We find that a Gaussian shape, though providing a reasonable first approximation, fails to account for the detailed shape of 

' M3's HB mass distribution. Indeed, the latter may have skewness and kurtosis that deviate slightly from a perfectly Gaussian solution. 

Alternatively, the excess of stars towards the wings of the distribution may also be accounted for in terms of a bimodal distribution in 
which both the low- and the high-mass modes are normal, the former being significantly wider than the latter. 
^ ' Conclusions. However, we also show that the inferred distribution of evolutionary times is inconsistent with theoretical expectations. 

This result is confirmed on the basis of three independent sets of HB models, suggesting that the latter underestimate the effects of 
\^ [ evolution away from the zero-age HB, and warning against considering our inferred mass distribution as definitive. 

I Key words, globular clusters: individual (M3) - stars: evolution - stars: horizontal branch 
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O ■ NGC 5272 (M3) is one of the best studied globiilar clusters in ^j^g^ ^^^^ ^3^^ predicted period distributions were most of- 

^ our galaxy. Extensive photoinetric studies of the cluster focusing ^^^^ ^^^^^^ ^^^^ observed. Among the possible solutions to 

on Its variable stars were published as early as the beginning of ^j^^ problem, these authors discuss the possibility of pulsation- 

the last century (Bailey 1913), and in spite of the many searches .^^^^^^ ^^^^ j^^^ ^^trapping" of HB stars at some 

^ . for ste lar variability in the cluster field that have been carried out .^^ evolution, as well as a rather peculiar multimodal 

since (e.g., Roberts & Sandage 1955; Szeidl 1973; Kaluzny et ^^^^ distribution (see §7 in Catelan 2004). 

al. 1998; Corwin & Carney 2001), until very recently previously „ . , , • ^ „ . ^ . „ 

unknown variable stars in the cluster field have continued to be ^ °" *^ "^^T"' Cartel am & 

reported in the literature (Clementini et al. 2004). Cassisi (2005) have recently argued, again on the basis of Monte 

Among M3's properties, of particular interest is its early Carlo simulations, that a sharply bmodal mass distribution may 

, f r>n T w indeed be able to account for the observed period distribution, 

classification, in terms of its RR Lyrae population, as an ^, , r , ■ , 

Oosterhoff (1939, 1944) type I (Ool) cluster. In fact, M3 is often ^he shapes of their two period distribution-based mass modes 

used to define the properties of the Ool class (e.g.. Smith 1995 turned out to be rather peculiar. More specifically their mode re- 

• f ■ \ r e 1 » » • 1 sponsible for the blue HB stars is essentially flat (between about 

and references therein). In spite of its role as the prototypical , ^ ^ , 

Ool cluster and its extensively studied properties, theoretical in- ^'^ V ° ^'^^ higher-mass mode responsible 

^ ^ ^- f ^, , ^- ^ ^ f\A^' DDT • ui tor the RR Lyrae and red HB Stars corresponds to a Gaussian dis- 

terpietation of the pulsation status of M3 s RR Lyrae variables ,, , ■ r^, , . 

■ ,u u- . f . T 1 D ^ ff tabution with <M) = 0.68 Mq and cr^ = 0.005 Mq. Thus their 

remains the sub ect of some controversy. In particular. Rood & , , , , . , , 

Crocker (1989) and Catelan (2004) have both called attention to suggested low- and high-mass modes are clearly detached from 

the fact that canonical stellar evolution/pulsation theory appears another Note, in addition, that Castellam et al. also advocate 

, . ui f f »u u 1 1 J J- * • truncating the high-mass Gaussian mode at its mean mass value 

somehow incapable of accounting for the sharply peaked distri- ? . , , , ^ . x . 

bution of (fundamentalized) periods that is found in the cluster '"^f ^"^"g *^ °f *^ Gaussian), in order 

(and possibly in other clusters as well). In particular, the Monte '° ^^"I'l overproduction of red HB stars. 

Carlo simulations presented by Catelan (2004), in which nor- As well known, the evolution of an HB star, for a given 

chemical composition, is primarily determined by its mass (or. 

Send offprint requests to: A. A. R. Valcarce more precisely, by the ratio between the envelope mass and the 
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core mass; e.g., Caloi, Castellani, & Tornambe 1978). Therefore, 
in the case of relatively simple populations as monometallic 
globular clusters (such as M3 itself; e.g., Sneden et al. 2004), 
one expects, within the canonical framework, the detailed HB 
stellar distribution on the color-magnitude diagram (CMD) to re- 
flect the underlying HB mass distribution. Accordingly, it should 
be possible, by careful examination of the available photometric 
data in such an extensively studied cluster as M3, to derive con- 
straints on the shape of the underlying HB mass distribution. 
This is of relevance in the present context, since - as akeady 
discussed - the HB unimodality that was adopted as a working 
hypothesis in the HB simulations computed by Rood & Crocker 
(1989) and Catelan (2004) has been identified as a possible ex- 
planation for the conflict between observed and predicted RR 
Lyrae period distributions in M3. In this sense, the main purpose 
of this paper is to check whether a multimodal mass distribution, 
as discussed by Catelan (2004), or a very peculiar bimodal mass 
distribution, as advanced by Castellani et al. (2005), is supported 
by the available CMD data. 

Unfortunately, not many studies have attempted to derive a 
mass distribution for M3 on the basis of the available CMD's. 
In this sense, while a normal distribution was favored by Rood 
& Crocker (1989), and while other authors have succeeded in 
reproducing several different HB morphology parameters of the 
cluster without the need to invoke a multimodal mass distribu- 
tion (e.g., Lee, Demarque, & Zinn 1990; Catelan, Ferraro, & 
Rood 2001b; Catelan 2004), data of much higher quality have 
become available in the literature since the late-80's for both 
variable and non-variable stars (e.g., Ferraro et al. 1997; Corwin 
& Carney 2001), thus warranting a more detailed look into the 
photometric evidence for bimodality (or lack thereof) along the 
M3 HB. This is the main goal of the present paper. 

In §2 we present the empirical data that were used in our 
study, along with the theoretical evolutionary tracks that were 
employed. In §3 we describe the method we used, along with 
tests that demonstrate that this method would be able to detect 
a bimodal mass distribution resembling the one suggested in the 
literature, if one were indeed present. In §4 we present the main 
results of our study. We close in §5 by presenting a summary and 
conclusions. 



2. Data 

2.1. Observational Data 

In order to reliably derive the mass distribution along the M3 
HB, we need a photometric database containing both variable 
and non-variable stars. The main problem we face in this enter- 
prise is related to the fact that, in order to derive reliable mean 
magnitudes and colors for the RR Lyrae stars, one needs time- 
series photometry covering the timespan of several days, and 
often weeks or months (particularly for variables with periods 
close to 0.5 d), whereas non-variable red and blue HB stars can 
be much more straightforwardly measured. Fortunately, and as 
already stated (§1), the M3 variables have been extensively stud- 
ied in the literature, and reliable quantities for their "equivalent 
static stars" have recently been provided by Cacciari, Corwin, & 
Carney (2005). 

For the variable stars, accordingly, we have used the Cacciari 
et al. (2005) database, which provides intensity-averaged mag- 
nitudes without the amplitude correction suggested by Bono, 
Caputo, & Stellingwerf (1995), since the magnitude of the equiv- 
alent static star, according to the latter's models, are always 
within 0.02 mag of the intensity-mean value (see Marconi et 



al. 2003). For the average colors, in turn, we have used the 
Cacciari et al. static colors, which do include the Bono et al. 
amplitude-dependent corrections. Accordingly, the total number 
of RR Lyrae stars in our study is 133, including 67 fundamental- 
mode (RRab or RRO) pulsators with regular light curves and 43 
presenting the Blazhko effect (RRbko), plus 23 pulsating in the 
first overtone (RRc or RRl stars). 

For nonvariable stars, in turn, we used the extensive photo- 
metric database by Ferraro et al. (1997), which includes BVI 
data for around 45,000 stars. For stars in the outer regions of 
M3 (r > 2'), they used both CCD data obtained at the 3.6m 
CFHT telescope and data obtained using photographic plates 
(Buonanno et al. 1994), whereas for the inner regions Hubble 
Space Telescope (HST) photometry obtained with WFPC2 was 
provided. In this work we have adopted their BV and VI data 
for M3's outermost and innermost regions, respectively. As dis- 
cussed in an Appendix, there may be a problem with the cal- 
ibration of their B-band (photographic) data, so that we have 
effectively utilized only their VI data to infer mass values for 
individual stars. 

2.2. Evolutionary Tracks 

We used the set of (canonical) HB evolutionary tracks from 
Catelan et al. (1998a) for a chemical composition Ims = 0.23 
(envelope helium abundance on the zero-age main sequence), 
Z = 0.001. These are the same evolutionary models that were 
used in Catelan (2004), and they also form the basis for the re- 
cent calibration of the RR Lyrae period-luminosity relation by 
Catelan, Pritzl, & Smith (2004). As a result, we have a total of 
22 evolutionary tracks for masses between 0.496 and 0.820 Mq, 
extending from the zero-age HB (ZAHB) to core helium exhaus- 
tion. 

The possibility of a spread in helium abundances, as re- 
cently suggested by several different authors in the case of 
very massive/peculiar clusters, such as NGC 2808, NGC 5139 
(w Centauri), NGC 6388, and NGC 6441 (e.g., Caloi & 
D'Antona 2007; Piotto et al. 2007), has not been taken into ac- 
count in this study, due to the lack of empirical evidence sup- 
porting such a hypothesis in the specific case of M3; in particu- 
lar, all the clusters for which a spread in Yus has been suggested 
present very long and well-populated blue HB "tails," whereas 
such a feature is lacking in the case of M3(3 In the present paper, 
we shall accordingly restrict ourselves to the canonical scenario. 

3. The Method 

The aforementioned evolutionary tracks are initially provided in 
the theoretical plane [log(L/Lo), log(reff)], whereas the empir- 
ical data are in the form of broadband filters-based magnitudes 
and colors. Therefore, in order to be able to infer the masses 
of HB stars on the basis of their positions in a CMD, we have 

' After this paper had been submitted, a new preprint surfaced sug- 
gesting that a (relatively small) helium excess may be present among 
the blue HB stars in M3 (Caloi & D'Antona 2008). Although not ex- 
plicitly noted by those authors, the color distribution along the HB 
is well known to be degenerate in terms of second parameter candi- 
dates (e.g.. Rood 1973), and so "vertical" information must be added 
in order to properly constrain the problem (e.g., Crocker, Rood, & 
O'Connell 1988; Catelan, Sweigart, & Borissova 1998b; Catelan 2005). 
Accordingly, in a forthcoming study we will apply extensive such tests, 
using both spectroscopy and photometric diagnostics, to quantitatively 
constrain the extent to which helium may be enhanced among the cooler 
blue HB stars of M3 (Catelan et al., in preparation). 
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Fig. 1. Evolutionary tracks for the adopted chemical composi- 
tion (solid lines), along with the corresponding ZAHB and he- 
lium exhaustion loci (dotted lines), are overplotted on the CMD 
for nonvariable stars (open circles) in the inner region of M3. 
The evolutionary tracks correspond to masses from 0.5 Mq (left) 
to 0.8 Mq (right), increasing in intervals of 0.025 Mq. The locus 
occupied by HB stars is schematically indicated by the dashed 
lines; stars within this region plotted with solid circles have had 
mass values inferred using our procedure. As can be seen, the 
adopted distance modulus, fiv - 15.00 mag, leads to a nice 
agreement between the observed and predicted lower envelopes 
of the HB distribution in this diagram. 



transformed the evolutionary tracks to the observational coor- 
dinates (Mb, My, Ml) on the basis of the color transformations 
and bolometric corrections by VandenBerg & Clem (2003) for a 
metallicity [Fe/H] ^ -1.5 and abundance of the alpha elements 
[a/Fe] ^ H-0.3 (e.g., Sneden et al. 2004). Interpolation on their 
tables was carried out using the algorithm by Hill (1982). 

In order to increase the internal precision in the mass deter- 
mination, we have generated additional evolutionary tracks, with 
a separation of 2 x 10^'^ Mq between consecutive values, by in- 
terpolating (also using the Hill 1982 algorithm) on the original 
tracks available to us. Note that Hill's is a very powerful Hermite 
interpolation algorithm, which indeed proved of great assistance 
in dealing with the non-linearities that are observed in the shapes 
of the evolutionary tracks as a function of mass. 

Apparent magnitudes of M3 HB stars have been trans- 
formed to absolute magnitudes using a distance modulus yuy - 
1 5 .00 mag in the V band. This is based on a comparison between 
our evolutionary tracks for Z - 0.001 and the M3 observations: 
as well known, HB stars spend most of their lifetimes close to 
the ZAHB, and this is only consistent with our theoretical mod- 
els for a distance modulus close to the indicated value. On the 
other hand, Harris (1996) provides jjy = 15.12 mag instead; if 
this value were adopted, we would find the rather untenable re- 
sult that more than 80% of M3's HB stars would be in a very 
advanced evolutionary stage. We later discuss the effect of un- 
certainties in fly upon our results. 
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Fig. 2. Evolutionary tracks for the adopted chemical composition 
(solid lines), along with the corresponding ZAHB (dotted line), 
are overplotted on the CMD for M3 RR Lyrae variable stars. 
Different symbols are used for the RRab (filled squares), RRc 
(empty squares), and RRbko (triangles) stars. The solid lines, 
from left to right, represent evolutionary tracks for 0.600, 0.625, 
and 0.650 Mq. A distance modulus jUv = 15.00 mag was again 
adopted. 



The colors were corrected for reddening using values for 
E(B-V) and E(V-I) of 0.010 (Harris 1996) and 0.016 mag (Rieke 
& Lebofsky 1985), respectively. The adopted E(B-V) value is 
only 0.003 mag smaller than the one implied by the Schlegel, 
Finkbeiner, & Davis (1998) maps. 

We inferred the mass and the evolutionary times (given as 
number fractions, where 0% corresponds to the ZAHB and 
100% to the helium exhaustion line ["terminal age HB," or 
TAHB]) for each individual HB star by choosing the track that 
most closely matched the star's observed CMD position (Figs.[T] 
and |2]l, with individual error bars being based on the observa- 
tional (photometric) errors (see Fig.[3]l. We have also carried out 
tests in which the evolutionary times are explicitly taken into 
account in the mass derivation 03.1.3I ). but have not found im- 
portant differences in the inferred mass distribution, at least in 
the case of M3-like HB morphologies, with respect to this more 
straightforward procedure. 

For the cases in which no individual error bars are provided 
in the original photometry, we have computed synthetic obser- 
vational error bars using the method described by Catelan et al. 
(2001b; see their §3.1), which invokes an exponential law for the 
errors in the observed magnitudes. The photometric error bars 
are of particular importance in the cases of stars lying close to 
the dashed line in Figure [1] since they define whether a star can 
reasonably be assigned to the HB phase or not (see also Fig. [3] 
right panel). 

In this sense, for stars falling at positions above the TAHB 
or below the ZAHB (dotted lines in Fig. 1), we used these obser- 
vational errors to determine, where applicable, the most likely 
mass value allowed for by the latter We considered that mag- 
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Fig. 3. Method used to determine the mass for each data point. Left panel: a red HB star The ZAHB is shown as a dotted line, 
whereas bold lines show evolutionary tracks for 0.65, 0.66, 0.67 and 0.68 (from left to right). The elHpse is the isoprobabihty 
contour for the 1-sigma error bars. The mass corresponding to the evolutionary track that passes closest to the center of the el- 
lipse {thin solid line, right) is 0.6744 M©, whereas that corresponding to the minimum evolutionary speed {thin solid line, left) is 
0.6726 Mq. The evolutionary lifetime fevoi along the HB is also shown in both cases, given as a percentage of the total HB lifetime. 
Right panel: a blue HB star Bold lines are evolutionary tracks for 0.545, 0.550, 0.555 and 0.560 Mq {from left to right). Here no HB 
track goes through the actual data point, and so we adopt for the star the mass corresponding to the track that comes closest to it, 
and still inside its 1-sigma error ellipse - in this case, 0.541 M©, which corresponds to a star at He exhaustion {thin solid line, left). 
The evolutionary track for the mass value with the smallest evolutionary speed that still goes through the 1-sigma ellipse is shown 
as a thin solid line (right); this corresponds to a star with a mass of 0.547 M© that has completed 99.8% of its HB evolution. 



nitude and color errors (cr^^ and (Xcoior, respectively) are stan- 
dard deviations of a normal distribution centered at the most 
probable value, which implies that a star has the same proba- 
bility to have a mass value determined by the points located at 
[My + (Tmv^ or [Mv, (B- V)o ± cTcoior]- In fact, one has 

concentric "isoprobabihty" ellipses centered at [My, {B-V)q\ 
(see Fig.[3]l. The mass assigned to a star thus corresponds to the 
track in our grid of evolutionary tracks intersecting the error el- 
lipse with the smallest semi-major axis. We used the ellipse with 
semi-major axis equal to the observational errors to assign er- 
rors to individual mass values {(Tm)- Note, however, that some 
stars can fall much above the TAHB or below the ZAHB lim- 
its, implying that they are not bona-fide HB stars within their as- 
signed photometric errors (i.e., at the 2 - cr level). For these stars, 
whose exact number depends on the distance modulus adopted 
(see §3.1.2 below), we do not attempt to assign individual mass 
values, since these would not be physically meaningful. 



3.1. Reliability Tests 

In order to ascertain the reliability of our adopted procedure, we 
have generated synthetic CMD's for the M3 HB, including syn- 
thetic photometric errors, with the purpose to verify whether our 
method is able to successfully recover a bimodal mass distribu- 
tion, if one is present, given the uncertainties in the photometry 
and in the distance modulus. Since for the synthetic models the 
mass distribution is known a priori, this provides us with a cru- 
cial test of whether the mass distribution to be inferred from the 



actual observations can be trusted, insofar as the presence (or 
lack) of bimodality is concerned. 

In a blind experiment, one of us (M.C.) computed the syn- 
thetic distribution and shifted the derived distributions using a 
distance modulus/reddening combination not known to the other 
(A.V.) - who in turn attempted to infer the corresponding mass 
distribution (see below). Reassuringly, the distance modulus and 
reddening favored by A.V. in the process were in excellent agree- 
ment (i.e., to within 0.01 mag in both magnitude and color) with 
the input values, which in turn were very similar to those given 
in the Harris (1996) catalog, namely: pv = 15.1, E{B-V) = 0.01. 

3.1.1. Synthetic Distributions 

We have constructed four input bimodal distributions using 
the Monte Carlo code sintdelphi (Catelan 2004 and references 
therein), as shown in the left panels of Figure|4] These mass dis- 
tributions are closely patterned after the Castellani et al. (2005) 
proposed bimodal mass distribution for M3 (see §1). Figure |4] 
{right panels) shows the recovered mass distributions, on the ba- 
sis of our method, assuming the same distance modulus as used 
when constructing the synthetic distributions (i.e., pv = 15.1). 
While it is clear that the added synthetic errors necessarily 
lead to some loss of information, and therefore to somewhat 
wider (inferred) mass distributions than in the synthetic mod- 
els, our method always succeeds in recovering a bimodal dis- 
tribution. Indeed, according to the KMM test (Ashman, Bird, 
& Zepf 1994), the distributions shown on the right-hand pan- 
els of Figure |4] are better described by bimodal rather than uni- 
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Fig. 4. Bimodality test: For four input bimodal mass distributions {left panels) we determined the mass distribution using our method 
(right panels). Due to the observational errors, we could not recover precisely the same input distributions, but in no case has a 
unimodal mass distribution been derived on the basis of an input bimodal distribution. 



modal Gaussian distributions, with probability always higher 
than 99.99%. 



3.1.2. Synthetic Distribution witii Variable Distance Modulus 

While the above results are suggestive, as we have seen there ap- 
pears to be some uncertainty in the M3 distance modulus, which 
may also impact the derived mass distribution. We have inves- 
tigated this source of systematic error by changing the adopted 
distance modulus over a wide range, from 15.0 mag to 15.4 mag, 
which should cover the full range of acceptable values (we recall 
that the value used in the simulations is 15.1 mag). In Figure|5] 



we show our results for three different input synthetic bimodal 
distributions (upper panels) for different values of the adopted 
distance modulus, from 15.0 mag (second row) to 15.4 mag 
(bottom row). As can clearly be seen, while the choice of dis- 
tance modulus does seem to affect the retrieved location of the 
peaks of the two output mass Gaussians, the bimodality in the 
mass distribution is always successfully recovered. Additional 
calculations and statistical tests using the KMM statistic con- 
firm that we only (mistakenly) retrieve a unimodal mass distri- 
bution when the adopted distance modulus is in error by more 
than about 0.3 mag. Note that, when a distance modulus value 
that is less than the correct one by over 0. 1 mag is assumed, we 



Valcarce & Catelan: Mass Distribution of HB Stars in M3 



-p-i— I— I— I— T-i— I— I— I— p-r 
INPUT 1 



200 

150 
100 
50 


100 - 
50 

100 
50 

100 
50 

100 
50 

100 
50 




I I I f I I I " I I I I I I I I 

Mv=l-5.0 
43.02% 




INPUT 2 



I I I 




INPUT 3 



Ai^=i5.o :: 

46.04% 



fj^= 15.0 
44.91% 




I I I 



(L4v=15.1 

100.00% 





f T r 



4=+ 




I -f I 1 ^ 1 I l-l 



H— f- 




(Uv=15.1 II 
99.43% 



T I — \ — h 

(Uy=15.2 
100.00% 




99.43% 




Mv=15.2 
100.00% 



jUv=15-2 
99.81% 



fP H- +11 



iLi„=15.3 
99.81% 




100.00% 



jUv=15.3 
99.81% 



+ I I 



-I— f- 



Aiv=15.4 
99.43% 



-I— f- 



99.81% 



_L 



Aiv=15.4 
99.43% 



_L 



0.6 0.7 0.8 

Mass [Mq] 



0.6 0.7 
Mass [Mg] 



0.8 



0.6 0.7 0.8 

Mass [Mq] 



Fig. 5. Distance modulus test: For three input bimodal mass distributions {upper panels) and assuming different distance moduli (as 
indicated in the insets, along with the percentage of recovered HB stars), we have derived the mass distribution using our method. 
Input distributions 1 to 3 are the same as in the previous figure. As can clearly be seen, we are always able to recover the input 
bimodahty, albeit not without incurring some systematic errors in the placement of the mass peaks. 



are unable to assign mass values to more than about 50% of the 
stars, because many will then fall below the ZAHB (see the last 
paragraph preceding §3.1). Naturally, such a loss of stars serves 
as an indicator that we have an incorrect distance modulus. In 
conclusion, our method seems to be very robust in its ability to 
detect mass bimodahty among HB stars in globular clusters. 



3.1 .3. Two Methods for the Determination of the IVIass 
Distribution 

While the method we used to infer masses basically looks for the 
evolutionary track that goes closest to the actual data point in the 
CMD, other methods can be devised in which the evolutionary 
times are used as a diagnostic criterion. In particular, to avoid an 
excess of stars close to the hehum exhaustion hne, where evo- 
lution is very fast, one can alternatively choose the evolutionary 
track that goes through the 1-sigma error elipse with the small- 
est evolutionary speed, thus in practice giving higher weight to 
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Fig. 6. Methodology test: For known input mass distributions (dashed lines) HB simulations with 1500 stars were computed as 
described in the text, and the mass distribution was then inferred [solid lines) from the resulting CMD distribution using two 
different methods: evolutionary tracks that pass closest to the individual data points (upper panels) and evolutionary track that 
passes through the 1-cr error ellipse with the smallest evolutionary speed Vevoi (bottom panels). On the left and right panels the cases 
of completely blue and completely red HB morphologies are shown, respectively, whereas the middle panel refers to the case of a 
more even distribution along the CMD, as more appropriate in the case of M3. As can clearly be seen, while the present method 
becomes less reliable for extreme HB types, it does provide an excellent description of the input mass distribution for an M3-like 
HB morphology. 



slower evolutionary stages. How would these different criteria 
affect the inferred mass distribution? 

In order to properly determine the evolutionary speed Vevoi, 
we first have to transform the CMD into a plane in which 
color and magnitude have comparable weights. This has been 
achieved using a technique similar to that described in Dixon 
et al. (1996), Catelan et al. (1998a), and Piotto et al. (1999), 
but here adapted to the My, V-I plane. Accordingly, we define 
rescaled "color" c and "brightness" b coordinates as follows: 



c = 168.3 (y-/)o + 96.64, (1) 

b = -42.67 Mv + 281.6. (2) 

With this definition, Vevoi is computed along an evolutionary 
track as 

VAc2 + Ab^ 

Vevoi = — ■ (3) 



Having thus defined the evolutionary speed, we have carried 
out numerical tests in which we generated, using sintdelphi, a 
synthetic CMD with 1500 synthetic stars and an input mass dis- 
tribution characterized by (M> = 0.642 Mq, tr^ - 0.020 Mq. 
Synthetic error bars were added as in Catelan et al. (2001b). We 
then checked the impact of the different criteria upon the inferred 
mass distribution. 

In Figure |6] we compare the mass distributions that were in- 
ferred using the two different criteria indicated in the beginning 
of this section with the input distribution. In the middle upper 
panel, the evolutionary track passing closest to the actual data 
point was selected. In the middle bottom panel, the evolution- 
ary track with the smallest evolutionary speed Vevoi within the 
1-sigma error ellipse was in turn selected. As can clearly be 
seen, the differences between the mass distributions derived us- 
ing either of these two methods are very small, and they are both 
clearly able to properly reproduce the input mass distribution. 
On the other hand, the reader is warned that the method does 
not perform as nicely for more extreme HB types. This is also 
revealed by Figure|6] where we show our attempts to recover the 
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Fig. 7. Mass distribution of nonvariable stars in the inner regions 
of M3. It displays two maxima, produced by stars pertaining to 
the red and blue HB domains. The minimum in between indi- 
cates the existence of stars in highly evolved stages towards the 
end of the HB phase. 




Fig. 8. Upper panel: inferred mass distribution for all RR Lyrae 
variable stars in M3. The mass distributions for each RR Lyrae 
subtype are shown in the other panels: RRab with well-behaved 
light curves {second panel), RRc (third panel), and RR Lyrae 
showing the Blazhko effect {bottom panel). 



input mass distribution in the cases of simulations computed for 
extremely blue left panel and red right panel HB morphologies. 
Irrespective of whether we adopt the purely geometrical criterion 
or the criterion involving the minimum Vevoi {upper and bottom 
panels, respectively), we are unable to reliably recover the in- 
put mass distribution in such cases. Photometry in other, more 
suitable bandpasses, where one finds a stronger dependence of 
colors and magnitudes on the stellar mass, would be required to 
achieve better results when the HB distribution is comprised of 
very blue or very red stars. 

We thus conclude that the simple method (in which the evo- 
lutionary track passing closest to the data point in the CMD is 
selected) appears to be good enough for our purposes, since it 
performs quite well for rather even, M3-like HB morphologies. 

4. Results 

4. 1 . Mass Distribution along tine HB of M3 

The mass distribution of stars along the HB of M3 was obtained 
separately for the nonvariable and variable stars, since we had 
to rely on empirical data in the VI bands for the former and in 
the BV bands for the latter. Our final M3 mass distribution corre- 
sponds to the sum of these two, separately derived, distributions. 

The non-variable stars that we selected for our study are 
those pertaining to the inner regions of M3. We inferred masses 
for 157 stars falling sufficiently close to HB evolutionary tracks, 
including 105 stars on the blue HB and 52 on the red HB {filled 
circles in Fig. [T]). The stars plotted with open circles within 
the dashed region of Figure [T] are those for which mass values 
were not assigned, since their photometry places them several cr 
away from the predicted HB locus. The resulting mass distribu- 
tion for these non-variable stars is bimodal, with peaks centered 
at 0.625 Ms and 0.665 Mq (Fig. [T]). The two peaks reflect the 



masses of blue and red HB stars, respectively; naturally, we an- 
ticipate that the gap in between these two mass modes will be at 
least partially filled when due account is taken of M3's variable 
stars. Had we used the method described by Rood & Crocker 
(1989), we would not have found any non-variable stars with 
masses around 0.640 Mq, since in that method the effects of evo- 
lution away from the ZAHB are not taken into account. 

The mass distribution for variable stars (Fig. [8] upper panel) 
was obtained from 127 of the 133 stars that we chose initially. 
The remainder of the stars fall below the ZAHB, and their pho- 
tometric errors do not allow them to be reconciled with the HB 
phase (some of them may be pre-ZAHB stars). A Gaussian fit 
provides a mean mass of 0.6445 Mq and dispersion 0.0076 Mq. 
More in detail, a Gaussian fit to the mass distribution for the 
RRab stars has mean mass and dispersion given by 0.6440 Mq 
and 0.0049 Mq, respectively (Fig. [8] second panel), whereas the 
Blazhko-type RR Lyrae have mean mass 0.6482 Mq and disper- 
sion 0.0085 Mq (Fig. [8] bottom panel). For the RRc, in turn, a 
rather uniform mass distribution is inferred instead, though in 
this case the number of stars is much lower (Fig. [8] third panel). 
Note, in any case, that the derived masses for the RRc stars tend 
to be lower than for the RRab stars, which is consistent with the 
expectations of canonical theory (see, e.g., Fig.|2]i. 

In order to study the global mass distribution along the M3 
HB, we must suitably combine our samples of variable and non- 
variable stars. In particular, we must ensure consistency with the 
observed proportions of stars falling along the blue HB, the red 
HB, and inside the instability strip. We accordingly use the pro- 
portions S : : = 39 : 40 : 21 (Catelan 2004), where S, 'V, 
and % indicate the numbers of blue, variable, and red HB stars, 
respectively. As a consequence, given the sample of red HB stars 
in our study, in order to derive an unbiased, final mass distribu- 
tion we must randomly remove 9 and 28 stars from our blue HB 
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and RR Lyrae samples, respectively. Our final derived mass dis- 
tribution is shown in Figure |9l 

We fit a Gaussian to this mass distribution, finding (Fig.|9]l 



(M) = 0.642 Mo, cr = 0.020 Mo 



(4) 



This result is consistent with that previously derived by Rood 
& Crocker (1989), who also find, using a method similar to 
ours but not taking into account evolutionary effects, a unimodal 
mass distribution, with (M) = 0.666 Mo and cr = 0.018 Mo. 
Catelan et al. (2001b) have also found that a unimodal mass 
distribution along the M3 HB is consistent with the observed 
HB morphology parameters of the cluster: according to their 
results, the innermost cluster regions may be characterized by 
(M) - 0.631 Mq and cr = 0.023 Mq, whereas the outermost 
regions, which seem to be a bit redder, may be described in- 
stead in terms of a normal distribution with (M) = 0.645 Mo 
ando- = 0.018 Mo. 

On the other hand, while the fit shown in Figure |9] seems 
rather acceptable, close inspection reveals what appears to be an 
excess of stars on the wings of the distribution, especially at its 
low-mass end. That the true mass distribution may be unimodal 
but not precisely Gaussian is also supported by an analysis of its 
skewness and kurtosis (see §14.1 in Press et al. 1992). We find 
values of Skew = 0.28 + 0.16 and Kurt = 3.31 ± 0.31, whereas a 
perfect Gaussian has Skew = and Kurt = 3. 

Naturally, since the derived unimodal distribution is not per- 
fectly Gaussian, an alternative description of the data can be 
accomplished with a linear combination of Gaussians. Indeed, 
when faced with the question to opt for a single Gaussian or a 
combination of two Gaussians, the KMM test gives, not surpris- 
ingly, preference to the latter, at a very high level of confidence 
(> 99.99%). The best-fitting bimodal solution is shown in Figure 
[TOl where the dashed line shows the result of the sum over the 
two derived Gaussians {solid lines). In this case, the best-fitting 
Gaussians are given by 



<M)i = 0.633 Mo, o-i = 0.026 Mo 
<M)2 = 0.650 Mo, cr2 = 0.008 Mo 



(5) 
(6) 



with 7 1 . 1 % of the stars in the first mode and 28.9% in the second. 

As previously discussed on the basis of Figures |4]and|5] one 
expects that the mass distribution determined in this way will dif- 
fer somewhat from the "intrinsic" mass distribution of the clus- 
ter. We have made an effort to account for this effect, at least 
in the case of the bimodal solution, by comparing the derived 
and input mass modes, as given by the panels with yUy - 15.1 in 
Figure|5] As a consequence, we infer that the "true" mass modes 
have on average cr values which are only about 77% of those 
given above. In like vein, the positions of the centers of the high- 
and low-mass modes are on average shifted by -0.008 Mo and 
-1-0.008 Mo, respectively. Therefore, the suggested (corrected) 
mass distribution for the M3 stars, assuming a bimodal solution, 
is the following: 



<M)i,,„, = 0.625 Mo, cri,coi- = 0.019 Mo, 

<M)2,cor = 0.658 Mo, (r2,cor = 0.006 Mo, 



(7) 
(8) 



again with about 7 1 % of the stars in the low-mass mode and 29% 
in the high-mass mode. 

Note that our derived mass distribution for M3 differs 
markedly from either the input or the derived solutions shown 
in Figures ID and |5] since the two Gaussians here are much closer 
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Fig. 9. Mass distribution for the M3 HB, obtained by combin- 
ing the mass distributions for the variable and non-variable stars. 
The best-fitting Gaussian (solid line) has (M) = 0.642 M© and 
a- = 0.020 Mo. 
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Fig. 10. As in Figure |9l but showing the best-fitting bimodal so- 
lution. The two Gaussians {solid lines) have (M)i - 0.633 Mq 
and cr, = 0.026 Mo (with 71.1% of the stars), and <M)2 = 
0.650 Mq and cri = 0.008 Mo (with 28.9% of the stars). 



together than was derived in almost all distributions except in 
those with the largest distance moduli (which were off from the 
correct solution by 0.2 dex or more in fiy). Therefore, while our 
study may give some support to the existence of some degree 
of HB bimodality in M3, the derived distribution also differs in 
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detail from the one that was set forth by Castellani et al. (2005). 
Interestingly, the two derived mass modes do not seem to be fully 
detached, contrary to what was suggested in the latter study. As 
a consequence, synthetic HB's based on the mass distribution 
that was inferred in this section are again unable to properly re- 
produce the observed fundamentalized period distribution in M3 
(see §1 for references to prior work on this topic). 

We repeated the whole procedure described in this sec- 
tion, but using the set of tracks independently computed by 
Pietrinferni et al. (2004). As a result, we derived a mass dis- 
tribution which is qualitatively very similar to that shown in 
Figures |9] and [TOl The only noteworthy difference is that, when 
using the Pietrinferni et al. tracks, the peak of the high-mass 
mode is slightly more pronounced than indicated in these plots. 

4.2. Evolutionary Times 

While we have previously argued 03.1.3I I that the inferred mass 
distribution is fairly robust in regard to the different recipes for 
using the evolutionary lifetimes to infer the mass values, we have 
so far not discussed what happens to the resulting distribution 
of evolutionary lifetimes itself. In principle, assuming a smooth 
feeding of the HB phase from the RGB tip, one should have 
an essentially flat distribution between and 100% of the HB 
lifetime, except of course for statistical fluctuations. How does 
the inferred distribution of evolutionary times look like, in the 
case ofM3? 

The answer is provided in Figure [TT](/e/0. In this plot, the 
upper panel shows the derived fgvoi distribution, corresponding 
to the mass distribution that was inferred by simply adopting the 
evolutionary tracks that go closest to each individual data point 
in the CMD. As can clearly be seen, there appears to be a dearth 
of stars in the range fgvoi ~ 20 - 60%, and/or an excess of stars 
with fevoi ~ 70-90%. The lower panel in the same figure reveals 
that, by adopting the smallest Vevoi value inside the 1-sigma error 
ellipse for each data point, not only is the problem not solved, but 
also a large excess of stars with fevoi ~ 10% results as well, as 
well as another strong peak of stars with fevoi ~ 55%. Figure [TT] 
{right) shows the result of the same exercise, performed using 
the Pietrinferni et al. (2004) evolutionary tracks. As can clearly 
be seen, the results are qualitatively very similar, thus indicating 
that the problem does not lie in our specific choice of evolution- 
ary tracks. 

As a matter of fact, the presence of these two peaks is 
an expected, direct consequence of the minimum Vevoi method. 
As well known, HB evolution is slowest close to the "turning 
points" on the CMD, including the blue nose and red "noses" 
that are clearly seen in the evolutionary track displayed in 
Figure|2l Therefore, any method based on minimum Vevoi values 
will tend to preferentially pick the fevoi values associated to these 
features - as in the case of Figure [3] (fe/f panel), where clearly 
the "blue nose" position was selected for a star whose original 
position on the CMD implied a slightly more advanced evolu- 
tionary stage. We conclude that the peak at fevoi ~ 55% can be 
ascribed to these "blue noses," whereas the peak at fevoi ~ 10% 
is instead due to the "red noses" that occur slightly after the star 
reaches the ZAHB. 

Irrespective of the method adopted, the main problem re- 
vealed by Figure [TT] is the fact that there are many fewer HB 
stars with 0% < fevoi < 50% than there are stars with 50% < 
fevoi ^ 100%. We have checked that this is not a problem af- 
fecting only a group of stars along M3's HB, but in fact the fevoi 
distribution looks bimodal for red HB, blue HB, and RR Lyrae 
stars. At present, we do not have an explanation for this problem. 
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Fig. 12. Luminosity evolution from the ZAHB, Alog(L/L0) = 
log L( fevoi) - log i( fevoi = 0), for three independent sets of evo- 
lutionary tracks. Evolutionary tracks whose ZAHB position lies 
close to the middle of the instability strip is shown in all cases. 
As can clearly be seen, according to the models one should ex- 
pect K. 65-70% of all stars along the horizontal part of the HB to 
lie within 0.05 mag of the ZAHB. Observations indicate instead 
that only » 40% are found within this range (see text). 



other than speculating that the present set of evolutionary tracks 
predicts too little luminosity evolution, thus leading to too few 
predicted stars at high luminosities, compared to the observa- 
tions (see also §2.2.4 in Catelan et al. 2001a, and §3.2 in Catelan 
et al. 2001b). We have checked that this is not a problem refer- 
ring exclusively to our adopted models; in fact, similar (or even 
more extreme) discrepancies are suggested by the Pietrinferni et 
al. (2004) or the Dotter et al. (2007) evolutionary tracks. This 
is further illustrated in Figure [12] which shows that, irrespective 
of the set of models used, HB stars are predicted to spend too 
little time at relatively high luminosities, compared with what 
is suggested by the observed CMD. More specifically, along the 
horizontal part of the HB, one expects to find, according to these 
models, » 65 - 70% of the stars within about 0.05 mag of the 
ZAHB. The observations, on the other hand, reveal ^ 40% of 
the HB stars within such a magnitude range from the ZAHB, 
even including in the tally those stars that fall below the ZAHB 
(presumably due to photometric errors; see, e.g.. Fig. |2]l. If we 
changed the adopted ZAHB position so as to better accomodate 
these "sub-ZAHB stars," the noted discrepancy would become 
even more dramatic. 

Naturally, until these problems are conclusively solved, mass 
distributions based on the CMD method, such as the one pro- 
vided in the present paper, should be considered tentative. In 
like vein, we caution that the fact that there are also too many 
(presumably) evolved red HB stars argues against a solution to 
this specific problem based on a component with a high helium 
abundance among the blue HB stars. 
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Fig. 11. Histograms showing the derived distributions of evolutionary times fevoi for stars in the HB phase, in terms of fractions of 
the total HB lifetime for the star's inferred mass (where 0% corresponds to the ZAHB and 100% to the TAHB). On the left we 
show the results obtained on the basis of our own evolutionary tracks, whereas on the right the results based on the independent set 
of tracks by Pietrinferni et al. (2004) are displayed. The upper panels show the evolutionary times based on the HB tracks that go 
closest to each individual data point in the CMD (e.g., those that would correspond to fevoi = 63.2% and 100% for left and right 
panels in Figure|3] respectively). The bottom panels correspond to the minimun evolutionary speed inside the 1-sigma error ellipse 
(e.g., those that would correspond to fevoi = 57.3% and 99.8% for the left and right panels in Figure[3] respectively). 



5. Summary and Conclusions 

In the present paper, we have studied the mass distribution of 
M3's HB stars by comparing their observed locations in the 
CMD with the predictions of evolutionary models. Our results 
suggest that, within the canonical framework, M3's HB mass 
distribution can be characterized either by a unimodal mass dis- 
tribution that is not perfectly Gaussian or by a bimodal mass dis- 
tribution in which the two mass modes are adequately described 
by Gaussians. In the latter case, the two mass modes appear to be 
separated in mean mass by only ^ 0.03 Mq, whereas the disper- 
sion in mass of the low-mass mode is significantly higher than 
that of the high-mass mode. 

As far as the distribution of evolutionary times is concerned, 
we find that it is not in good agreement with the canonical ex- 
pectations, in that we obtain a bimodal distribution, with a dearth 
of stars with fgvoi ~ 20 - 60% and/or an excess of stars with 
fevoi ~ 70 - 90%. This suggests that the present evolutionary 
models underestimate the luminosity evolution along the HB. 
We have checked that other sets of evolutionary tracks, such 
as the ones computed by Pietrinferni et al. (2004) or Dotter 
et al. (2007), indicate similar, or even smaller, luminosity evo- 
lution than the present ones. Until this problem is solved and 
better agreement between observed and predicted fevoi values 
can be achieved, mass distributions derived using a CMD-based 
method, such as the one presented in this paper, should be con- 
sidered tentative. Conversely, no solution for the HB star prob- 
lems in M3 that were previously discussed by several authors 
(Rood & Crocker 1989; Catelan 2004; Castellani et al. 2005; 
D' Antona & Caloi 2008) can be considered complete until good 
agreement between predicted and observed lifetime distributions 
is finally achieved (see also Catelan et al. 2001a,b). 
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Appendix A: Photometry in the BV Bands: Is there 
a Problem? 

In Figure lATI we overplot our evolutionary tracks with the obser- 
vational data for the outer regions of M3. As can be seen, there is 
a clear disagreement between the predicted locus of M3's blue 
HB stars and the observations. We find no combination of fiy 
and E(B-V) that allows us to provide good agreement between 
the models and the empirical data. On the other hand, for the in- 
ner regions of M3 (Fig.[TJ no such problem was present. Since 
the measurements for the inner regions are based on CCD data 
(Ferraro et al. 1997) whereas for the outer regions photographic 
data were used (Buonanno et al. 1994), this strongly suggests 
an error in the calibration of the photographic data, especially in 
regard to the color term. This is consistent with the discussion 
in Ferraro et al., who have also called attention to the possibility 
of large errors in the B-band photometry. Note that this problem 
does not affect the adopted colors and magnitudes for the RR 
Lyrae stars, whose adopted BV photometry (§2) came from an 
independent source (namely, Cacciari et al. 2005 and references 
therein). 
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Fig. A.l. My, B-V CMD of nonvariable stars in the outer regions 
of M3. Evolutionary tracks (solid lines), along with the corre- 
sponding ZAHB and TAHB lines (dotted lines), are overplotted. 
As can clearly be seen, for colors bluer than (B-V)o = 0.2 the 
models progressively deviate from the empirical data. As dis- 
cussed in the main text, this is likely due to an error in the cal- 
ibration of the photographic data [compare with Fig. 1, which 
reveals excellent agreement between the models and the CCD 
data in the My, (V-I)o plane]. 



